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Abstract 

We calculate smoothed correlators for a large random matrix model 
with a potential containing products of two traces tiW\{M) -trW^M) 
in addition to a single trace trV(M). Connected correlation function 
of density eigenvalues receives corrections besides the universal part 
derived by Brezin and Zee and it is no longer universal in a strong 
sense. 
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1 Introduction 



Random matrix theory was originally introduced by Wigner and studied in 
detail by Dyson and Metha to investigate statistical properties of energy 
levels of heavy nuclei QXJ] . It has been applied to various fields recently such 
as quantum chaos ||, quantum dots or 2d discretized gravity ||. It has also 
been applied to quantum transport problem of mesoscopic wires reviewed in 

I- 

In the above mentioned contexts universal behavior of correlation func- 
tions of eigenvalues has been discussed. There are two types of universali- 
ties, one for a short distance behavior (where correlation functions oscillate 
rapidly) |J and the other for a smoothed correlator of distance scale larger 
than the rapid oscillation ||. In this letter we study the latter universal 
behavior of correlation functions. 

The universal large scale behavior of density correlation functions was 
pointed out by Brezin and Zee who stressed its importance in disordered 
systems though it was already derived by Ambjorn, Jurkiewicz and Makeenko 
[0] in the context of 2d discretized gravity. They considered a random matrix 
theory for large matrices M, whose statistical weight is given by 

}_ e -NtrV(M) /-Q 

and showed that, in the large matrix limit, the connected two point correla- 
tion function of density of eigenvalues (p(x)p(y)) c has a universal form which 
has no explicit dependence on the potential V. For a symmetric potential it 
is 

1 a 2 — xy 



2n 2 X(x - y) 2 ^( fl 2 - x 2 )(a 2 - y 2 ) 



(2) 



where a and —a are the end points of the density distribution and A = 
1/2, 1, 2 corresponding to orthogonal, unitary or symplectic ensembles. This 
universal form has been calculated in various methods |7|, [6|, §, |9|, [12]. and 
its 1/N corrections were also studied for general ensembles [7|, |10 , 



Physical implication of this universality to the universal conductance fluc- 
tuation of mesoscopic wire is discussed in a recent review by Beenaker [|]]. 
Universality here means that fluctuation of conductance 5G in a metalic 
regime is of order e 2 /h, independent of sample size or disorder strength. 
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This strong suppression of conductance fluctuation was first pointed out by 



Altshuler 113] and Lee and Stone [14|. Also, when magnetic field is applied, 



variance (5G) 2 of the conductance decreases exactly by a factor of two com- 
pared to the fluctuation without magnetic field. The system can be modeled 
by random matrix theory for transmission matrix ||]. Since the variance 
of conductance is proportional to the connected density-density correlation 
function, the above universality of conductance fluctuation is mapped to the 
universality of correlation functions in random matrix theory. Looking back 
to the equation @, first it is 0(N°) and independent of the system size (N) 
or details of disorder (V). Next it decreases exactly by a factor of two when 
magnetic field is applied (A is changed from 1/2 to 1 by applying magnetic 
field.) It is quite nice that such a simple model as the random matrix theory 
can explain some of important features of complicated systems. 

Since the main guiding principle of random matrix theory is randomness 
and symmetry, it is natural to ask whether this universality still holds for 
more general ensemble which is invariant under symmetry rotation. Simple 
generalization of the potential ([[]) is to add products of traces to the statistical 
weight of matrices; 

}_ e -NtvV(M)-tvW 1 (MytrW 2 (M) /g\ 

Z 

This ensemble is invariant under M — > \JMU~ 1 where U is an orthogonal, 
unitary or symplectic matrix correspondingly. The ensemble with this gener- 
alized potential was studied by fl5j in the context of 2d gravity. Universality 
in this ensemble was discussed by [16, I7| . Brezin and Zee [HJ argued that 



this model is equivalently described by an ensemble with an effective single 
trace potential V e ff and concluded that the universality still holds. 

In this letter we study this generalized ensemble and obtain density of 
eigenvalues and its correlation functions explicitly. In section 2, we review a 
collective field theory approach to an ensemble with a single trace potential 
and show how the universal behavior of correlation functions emerges. In 
section 3, we generalize it to an ensemble with a multi trace potential (con- 
taining products of two traces). We show that the universality is broken and 
the correlation function is no longer universal in the strong sense. Finally in 
discussion, we discuss why the argument by Brezin and Zee [16] does not hold 
for correlation functions. In Appendix A and B, we prove useful formulas 
and in Appendix C we give an example of a correlation function for a simple 
ensemble with a multi trace potential. 



2 



2 Single Trace Matrix Model 



In this section we review how to calculate density of eigenvalues and its 
correlations in random matrix theory in the large N limit (N is a size of 
matrices) and show how the universal form of a two-point correlation function 
emerges. We consider a matrix model with an ordinary single trace potential 
in the Collective Field approach [Tj|. This approach is easily generalized to 



a multi trace potential, discussed in the next section. 
The free energy F[V] is defined as follows: 

e -Ni F[ V] ^l d J± e -NtrV(M) = f dxi _ dxN A 2A V(*0 (4) 

J Vol J 

where Vol is the volume of gauge symmetry group, A = ni<i<i<jv( x i ~ x j) * s 
the Van der Monde determinant, and A = l/2,lor2 for orthogonal, unitary 
or symplectic ensembles correspondingly. Partition function is invariant un- 
der orthogonal, unitary or symplectic rotations and the matrix integral can 
be reduced to integrals over its eigenvalues x± (i — 1 ~ N). (Normalized) 
density of eigenvalues is defined by 

1 N 

P( X ) = wJ2 S ( x - x i)- ( 5 ) 

iv i=l 

Connected density-density correlation functions can be obtained from the 
free energy by taking variational derivatives with respect to V(x): 

It is obvious now that the leading non- vanishing term for an m-point (unnor- 
malized) density correlation function is 0(N 2 ~ m ). The standard procedure 
of collective field theory is to rewrite integrals over eigenvalues in terms of a 
functional integral over density p{x). Inserting one to the equation ([|) 

1 = jDp(x)l[6^p(x)-^j25(x-x^, (7) 

we obtain up to an overall constant (see Appendix A) 

e - N 2 F [ V ] = J Dp ^ da e N 2 X J dxdyp(x)ln\x-y\p{y)-N 2 J dxp(x)V(x)+ia( J dxp(x)-l) 

(8) 
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Lagrange multiplier a is introduced to impose the constraint that an integral 
of p(x) is normalized to one. We are interested in the large N behavior and 
the Jacobian J[p] can be neglected in this limit (see ]nj for details). The 



resulting integral over p can be evaluated by steepest descent method which 
requires to solve the following equations of motions ]18| : 



= 2N 2 \Jd£ln\x-t\po(0-N 2 V(x)+ia , (9) 
= J dx p {x) - 1. (10) 
By differentiating the first equation we have the equation of BIPZ [21 

(id 

and it determines the stationary value po(x). We assume here that po(x) is 
equal to zero for x < a or x > b (one-cut from a to b). Then the solution of 
this Cauchy integral equation is given by |20[| : 



Po(x) = - p L dt ^ - *V(fl + 1 (12) 



n J(x-a)(b-x) \ 2rrX ' a € 

The second term in the curly bracket is a solution of a homogeneous integral 
equation (i.e set the r.h.s. of eq. (|LTD zero) and its coefficient is determined 
such that po(x) satisfies equation (|I0| ), while the first term does not contribute 
to the integral of eq. (|T0"D (see Appendix B). 

Equation (0) defines cr [V] as a functional of V and, of course, independent 
of x even though it enters to the solution manifestly. Unknowns left are the 
positions of the end points of the cut [a, b]. By choosing boundary conditions 
Po(a) = and po(b) = we get: 



p„(a) = _!_/ ^/LJL/te) + 1 = 0, (13) 




Po (6) = ^ -^. id ^izj V ' i0 +1= „. (14) 

Equations (|l3l) and fll4|) determine a = a[V] and 6 = 6[V] as functionals of 
V(x), and should be solved first. This completes the solution po(x) through 
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the equation (0). For a simple case V(x) = x 2 /2 this procedure reproduces 
the famous Wigner semi-circle solution. 

Dependence of Pq{x) on V(x) comes both explicitly, and implicitly through 
the end points and in order to find connected correlation functions @ we 
have to know V^-dependence of the end points. However we have an impor- 
tant relation that, given the boundary conditions ( [13|) and ([14]), the solution 



(|12 ) satisfies 

d d 

d~a Mx) = db Mx) = °' (15) 
This can be checked by straightforward calculations. Therefore upon varia- 
tion of po(x) with respect to V(x) only manifest dependence on V(x) in (|T^) 
is relevant, while V^-dependences of the boundaries a and b are cancelled out. 
Since the two-point correlation function is obtained from po(x) by taking a 
variational derivative, this is essentially statement of the universality of a 
two-point correlation function: 



& 1 d 

N 2 {po(x)p (y))c = - 7777 , po(x 



//( (y-a)(b-y) P 

SV(y) 2tc 2 X dy\\ (x — a) (b — x) y — x 



(16) 



All dependence on the potential V(x) is implicit only through the end-points 
a and b. Here P stands for the principal value: P^ = \im e ^ox/(x 2 + e 2 ). 
This result has been obtained in various papers by various methods || [7]. ^|. 

IBB 



Collecting all results together we obtain for the free energy F[V] the 
following expression (in the N — > 00 limit): 

fb[V] fb[V] 

F\y] = / dx po(x)V(x) — A / dx dypo(x)ln\x — y\ po(y). (17) 

Ja[V] Ja[V] 

As expected, density of eigenvalues is given by the saddle point solution; 

SF[V] 



(p(x))< 



SV{x) 



po(x) + J d£ J^j(V(0 - 2A / <hj h,\.r >)\ PnHj)) 



= Po(x) - ^ J_ d£ (po(Opo(x))c = Po{x). (18) 

In principle 1/N corrections can be calculated by evaluating the Jacobian 
J[p] [[LJ| and fluctuation around the saddle point. 



b[V] 
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3 Multi Trace Matrix Model 

Let us now generalize our approach to a multi trace case. We consider an 
ensemble with the following statistical weight 

iexp ( - NtiV(M) - ~ £ u a ptrW a (M)trW p (M)). (19) 

Z V 1 aJ3=l J 

In this case, after integrating over the angular variables, free energy is given 
by: 

e -N*F[V] = J dxi _ dxN A 2A e -iVEtl^)-|Ef l/3=1 ^E5=i^^)^fe) 

(20) 

and as before connected m-point correlation functions can be obtained by 
taking functional derivatives with respect to V{x) (see (^)). Expressing F[V] 
in terms of the density p(x) we obtain: 

e -N 2 F[V] = J JJp^dfj e N 2 \ J dxdyp(x)ln\x-y\p(y)-N 2 f dxp{x)V(x) x 
a, l 3=i UJ °'l3 f dx p(x)W a (x) J dy p(y)Wp(y)+ia( J dxp(x) l) j^j 

with the same Jacobian J[p] as in the previous section. Again in the leading 
order in N we can set J[p] = 1. The steepest descent equations are: 

K 

= 2N 2 X / d£ln\x - f |p (0 - iW(ar) - N 2 ^apW a (x)c p + ia , 

a,p=l 

(22) 

= J dx p {x) - 1 (23) 
where cg's are constants 

op = J dx p (x)Wp(x) (24) 

which are determined self-consistently later. Taking a differentiation of the 
equation (22) we have a generalized equation of BIPZ 

p /«£S=^ + s?, (25) 



a, /3=1 
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and it determines the stationary value po(x). Let us consider again a one-cut 
solution, i.e po(x) is equal to zero if x < a or x > b. Then the general 
solution is given by |2"0[] : 

Po(x) = [ b d{G(x,Z){v(Z) + £) o>apW a (£)cf\+ - L = (26) 

7a 1 «,/?=! J *yj(x-a)(b-x) 



where G(x,y) is a differential operator defined by 

f 



G(x,y) 



(y-a)(b-y) P \ d_ 
27T 2 X I \ {x — a)(b — x) y — x J dy 



(27) 



The coefficient of the homogeneous part (the second term) of po(x) is deter- 
mined so as to satisfy the constraint fl2"3"|). 

The constants cg's are still unknown. In order to fix them we plug equa- 
tion (|26| ) to the equation ( p^) and obtain a set of algebraic equations for cp's, 
which can be written in the following compact form: 



A 



where 



O 



7=1 



6 „ 

dxdyWp(x)G(x, y)V(y) + / <ir- 



a 7r^/ (a; — a) (6 — x) 
fy^ = 8p 1 -^u ai j dxdy Wp(x)G(x,y)W a (y). 

a=l Ja 



Assuming further that 



det |0| ^ 



equation (2^) can be inverted 

c a = Y^i^^^Op 

to give the solution for po(x): 

K 

Po[x) " " ' ' ' 



(28) 

(29) 
(30) 

(31) 
(32) 



/ d£G(x,o\v(0 + £ W a {i)uo a ^- l ) M 



+ 



7T 



x — a)(b — x) 



(33) 



Now the only unknowns left are the end points of the cut [a, b\. We can 
fix them by choosing boundary conditions po(a) = and po(b) = 0. These 
equations determine a = a[V, W] and b = b[V, W] as functionals of V(x) and 
W(xYs. Then equation (]33|) provides final expression for po(x). 

As in the previous section one can prove that the solution fl33|) for po(x) 
satisfies 

^p (x) = 0. (34) 



In order to prove it, notice that, from equation ( j26|) and the boundary con- 
ditions po(a) — po(b) = 0, we get 

^ P0 (X)= fdZG(x,0 J2 WofiWaiO^Cfi. (35) 
a a, 3 1 

All other contributions cancel out due to the boundary conditions at the 
end points. Therefore taking derivative of cp in equation ([24|) we get a set 
of homogeneous algebraic equations, which can be written in the following 
form 

E^^ = 0- (36) 
Since det \Q\ ^ 0, as previously assumed in equation (|31|), we conclude that 

w = ° (37) 

and equation (|3^) is proved. Similarly one can prove that 

dp (y) dc a 



da da 



0. (38) 



Therefore upon variation of po{x) with respect to V(x) only the manifest 
dependence on V(x) in ( |3"3"D remains. For a connected two point correlation 
function we obtain: 



AT 2, ( W u 5 Po(y) 1 ^ f {y-a){b-y) P 
N {p{x)p{y)} c = 



SV(x) 2tt 2 X dy\\ (x — a)(b — x) y — . 
J2 f 'dZG(x,0W a (0 a aP [ b dCG(y,C)W f5 (0 (39) 

a 1 ■/ a Jo 



a,/3=l 
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where 



A' 



er, 



£ u; aP (Q- 1 )p T (40) 



p=i 

Correlation function Q53) is symmetric for x and y which follows from the 



symmetricity of the matrix u and the definition of f2 (30|). 

The correlation function (|39"D manifestly depends on W a 's and no longer 
universal in a strong sense. However it is independent of V(x) and also the 
short distance behavior is dominated by the first term of the universal form, 
which is what we call weak form of universality. 3 It is 0(N°) and there is still 
strong suppression of fluctuation though its amplitude is not universal. And 
the second term of the correlation function which breaks the universality 
is not inverse-proportional to A generally. But if the second term in the 
definition of Q eq.(|50"|) is much larger than one the correlation function ( |59]) 
is again inverse-proportional to A. 

The free energy F[V] is written in the following form: 

r b[V,W] K 

F[V] = / dx V(OPo(0 + E ^afiCaCp 

Ja[V,W] a> p- x 

- A / dxdyp (x)ln\x - y\p Q (y) (41) 



where c a 's are defined in (p2|) . 

Finally we consider that LUi2 = UJ21 = 1 and all other components 

equal to zero. Using equation fl30D and (140|) the explicit form of the two-point 
correlation function is 



N 2 (p(x)p(y)) c ' ° 



(y-a)(b-y) P 



2n 2 \ dy\\ (x — a)(b — x) y — x . 
? dtG{x,t)W x {£) {W2 '^ 2) I" dCG(y,C)W 1 (0 

J a det J a 



3 In the paper by B. Eynard and C. Kristjansen |2^| they proved universality of cor- 
relation functions for O(N) model on random lattice which can be written as one matrix 
model with an infinite sum of products of two traces. The universality might be recovered 
for their case since the infinite sum may make the determinant eq. ( |3l| ) divergent and a aj 
in eq.(p0|) zero. We would like to thank C. Kristjansen for calling our attention to their 
papers. 
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h - i _ (Wi ■ G ■ Wo) r b 

d^G(x,OW l (0 *± j a d(G(y,()W 2 (0 

b i- (Wo ■ G ■ Wt) r b 

d£G(x,OW 2 (0 Wl) j a dCG{y,C,)W x {0 

d£G(x,OW 2 (0 { 1 lJ / d(G(yX)W 2 ((), (42) 

det Ja 

where 

det |0| = (1 — (W 2 ■ G ■ W^il — (Wi ■ G ■ W 2 )) - (W x ■ G ■ W 1 )(W 2 ■ G ■ W 2 ). 
(W a -G- W p ) = fdxdy W a {x)G{x,y)W p {y). (43) 

J a 

A concrete example for W\ = W 2 = g 2 x 2 /2 + (? 4 x 4 /4 is given in Appendix C. 

4 Discussion 

In this letter we studied density-density correlation functions for random 
matrix models with a generalized potential containing products of two traces 
trWi(M) • trW 2 (M) in addition to a single trace trV(M). We showed that 
the two point function has no longer the universal form which depends on 
the potentials only through the end points. This is against the argument by 
Brezin and Zee [0. They have considered a special ensemble 

p(M) = ! e -^trV(M)-[W(M)] 2 /2_ ^ 

They argued that the effect of the second term is just to renormalize the 
potential V to V + a W where a is a constant determined by V and W and 
claimed that the universality for two point correlation function still holds for 
the above ensemble. Here we briefly review their argument and discuss why it 
cannot generally hold for higher point correlation functions. By introducing 
an auxiliary variable a the partition function is written as 

Z(V,W) = e - N2F[v ' w] = J dMe- {NtrV{M)+[tTW(M)]2/2) 

= J da n^nv^ (45) 
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up to an irrelevant overall factor and the integral over a runs over the imag- 
inary axis. In the large N limit, the integral over a can be evaluated at the 
saddle point and 

2 

F[V, W) = F[V + a W, 0] - ^ (46) 
where ao is determined by 

a = ^-F[V + aW J 0}\ ao . (47) 

ato depends on the details of the potentials. We can now obtain its density 
distribution function: 



SF[V, W) 
5V(x) 

SF[V + a W, 0] . ,dF[V + a W,0] . 5a 
+ ( aZ. «o)- 



6V(x) ]a ° ' v da ~ UJ 5V(x) 
5F[V + a W, 0] 



|a - 



(48) 



5V(x) 

This density distribution is the same as that of an ensemble with an effective 



potential tiV e ff = tr(V + a W) as discussed in jig] . In order to obtain two 
point correlation function we then take variational derivative of the distribu- 
tion function with respect to the potential V(y). po(x) depends on V(y) not 
only explicitly but implicitly through a . It becomes 

5 2 F[V,W] 5 Po (x) 



G(x,y) 



5V(x)5V(y) 5V(y) 

S 2 F[V + a W, 0] dp (x) 5a 



5V(x)5V(y) 10 da 5V(y)' 



(49) 



The first term is the universal correlation function for an effective potential 
V e ff and actually independent of the details of the potential. But the second 
term, which is a product of a function of x and that of y, depends on the 
potential explicitly and the universality is broken down. (It is straightforward 
to show that the above expression is equal to a special case of that obtained 
in section 3. ) In this simple case of potential, the extra term is factorized 
into functions of x and y and we may say that there is still universality in 
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a weak sense. As we discussed in the paper, if the potential term is more 
complicated, other products of functions at x and y are added and this weak 
universality is also gradually broken. 

Our result can be also generalized to potentials containing products of 
more than two traces. Two point correlation function is again given by the 
universal form plus a sum of products of a function of x and that of y. In 
this case we have to solve a non-linear equation to obtain these functions. 

To conclude we have shown that the two point correlation function is 
no longer universal in the strong sense and it depends on the details of the 
potentials. This implies that, when we apply it to the random matrix models 
for transmission matrix of a quantum wire, conductance fluctuation is not 
exactly universal. Amplitude of the fluctuation might depend on the system 
size or disorder strength. Also the variance of the conductance does not 
decrease exactly by a factor of two when magnetic field is applied. 
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Appendix A 

We prove the following identity we have used in equation (H): 

N r+oo 

ln(A 2 ) = ln\xj — Xj\ = N 2 / dxdy p{x)ln\x — y\p(y) + const. (A.l) 

A-LA J —CO 

Regularizing 5-function in the definition of the collective coordinate p(x) (see 
(15D) as follows, 



p(x) = lim ±= £ -. ^ ( A -2) 
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we obtain 

/+oo 
dxdyp(x)ln\x — y\p(y) = 
-oo 

= hm> / cbcfy- Zna-y A.3 



'■■J 



The sum can be separated into two pieces: i ^ j and i — j. The % ^ j terms 
lead to the expression for the Van der Monde determinant. The i = j terms 
become, by setting £ = x — Xj and ( = y — X{, 

e / 7r e / 7r 



N f+°° e/n 

limV / dxdy 7 - j/ 

e ^~{J-oo (x — Xi) z + e z 



x — Xj) 2 + e 2 (y — Xj) 2 + e 2 

= £&E HZ - CI 7^ = const (A.4) 

and equation (|A.1| ) is proved. 



Appendix B 

Due to the following identity 

f(y) = - f dx 1 =-^- = 0, a<y<b, (B.l) 

nJa J(x-a){b-x)V- x 

the inhomogeneous term in po(x) (equation (|12|) or equation (|26|) ) does not 
contribute to an integral (|10D. This can be proved as follows. We first define 
a function 

1 r b 1 1 

F(z) = - / rfx . = . (B.2) 

nJa J( x -a)(b-x)z-x 



Then it follows 



f(y) = l[F(y + ie) +F(y-ie)), e - 0. (B.3) 



Choosing the square root to be positive on the upper side of the cut (and 
negative on the lower), the line integral ( |B.2| ) can be deformed to a contour 
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integral along a path C circling clockwise the cut between a and b. z is 
outside of this contour. Since the integrand vanishes at infinity the contour 
can be deformed smoothly to wind counterclockwise around z; 



F(z) 









j> dx 












j> dx- 


2^ 





1 



(x — a)(b — x) z x 
1 1 

x — a) (b — x) z ~~ x 



y/{z-a)(b-z) 



(B.4) 



F(z) has different signs on the upper and lower sides of the cut and we get 
f(x) = 0. 



Appendix C 

In this appendix we give an example for an ensemble with a potential NtrV(M) + 
(trW(M)) 2 where 



W{x) 



g 2 x 2 g A x A 



(C.l) 



2 4 

We assume that V(x) is a symmetric potential and a = —b. From eq. (ff2]) 
we have 



N 2 (p(x)p(y)) c 



1 d f b 2 -y 2 P 



2ir 2 \ dyY\b 2 — x 2 y — x 

> 6 ^«^« ,-2 { w-G.w) L d<6(y ' 0WKy (c ' 2) 



First it is straightforward to evaluate the integral 



" d£G(x,0W(0 = 2 , L g/ ^ 

-6 2-K 1 \\Jb 2 - x 2 J-b 



w'(0 



2rr 2 XVb : 
1 



[ b dtdb 2 - £ 2 
'b 2 - x 2 J-b ^ V ^ 



P P 

g 2 (l + x- — -) + g A (£ 2 + x£ + x 2 + x 3 



£ — x' 



2<!r 2 \y/b 2 - x 2 



Tib 2 



92{— ~ ttx ) +g±( 



,7r6 4 irb 2 x 2 



TTX 



£ — X 

(C.3) 
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Then by using 



2m 



(2m - 1)!! 
2mH 



nb 



2m 



we can obtain 



(W-G-W) 



9%b A , g 2 g±b 6 i 9 268 



32 



+ 



32 



+ 



1024' 



(C.4) 



(C.5) 



The end points —b and b are determined by the conditions that the density 
distribution should vanish at end points. The final answer for the connected 
density-density correlation is in general not inversely proportional to A. This 
means that, when we apply matrix theory to the problem of conductance 
fluctuation, the variance (5G) 2 does not decrease exactly by a factor of two. 
However, if \(W ■ G ■ W)\ » 1, we can approximately write the two point 
correlation function as 



N 2 (p(x)p(y)) c 
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2n 2 X dy IV b 2 — x 2 y — x 
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ny ) 
(C.6) 



and it is again inversely proportional to A. 
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